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Topological phase space disconnection has been recently found to be a general phenomenon in 
isolated anisotropic spin systems. It sets a general framework to understand the emergence of ferro- 
magnetism in finite magnetic systems starting from microscopic models without phenomenological 
on-site barriers. Here we study its relevance for finite systems with long range interacting potential 
in contact with a thermal bath. We show that, even in this case, the induced magnetic reversal 
time is exponentially large in the number of spins, thus determining stable (to any experimental 
observation time) ferromagnetic behavior. Moreover, the explicit temperature dependence of the 
magnetic reversal time obtained from the microcanonical results, is found to be in good agreement 
with numerical simulations. Also, a simple and suggestive expression, indicating the Topological 
Energy Threshold at which the disconnection occurs, as a real energy barrier for many body systems, 
is obtained analytically for low temperature. 

PACS numbers: 05.20.-y,05.10.-a, 75.10.Hk, 75.60.Jk 



I. INTRODUCTION 

One dimensional toy models are widely investigated 
in statistical mechanics 1]. From one hand, the possibil- 
ity to get analytical results represents the starting point 
for analyzing more physical models. On the other hand, 
due to their high simplicity, they allow a better under- 
standing of the key mechanisms at the basis of impor- 
tant physical effects. It is the case of the Topological 
Non-connectivity Threshold (TNT), recently introduced 
and addressed in f3] and investigated in other related pa- 
pers 3, 4, 5]. In these simple toy models, with a well 
defined classical limit two key features were introduced, 
anisotropy and long-range coupling. Even if acting in dif- 
ferent ways, they are both essential to generate a signifi- 
cant disconnection of the Hamiltonian phase space lead- 
ing to what is known in literature as breaking of ergod- 
icity d, @ . While anisotropy is a common paradigm in 
the phenomenological models of ferromagnetism (usually 
introduced as on-site anisotropy barrier in microscopic 
models) Q, long range interactions were re-discovered 
quite recently, due to the development of powerful and 
efficient techniques 8]. Indeed, strictly speaking, a well 
known model for anisotropy including a — term in 
the magnetic energy [9] (M ^ being the sum of 

all magnetic moments within a suitable domain) exactly 
matches an all-to-all interacting model close to what we 
consider here below. The role of anisotropy in finite spin 
systems has attracted much attention recently, following 
the experimental verification of ferromagnetic behavior 
in finite ID systems with strong anisotropy [TJ], con- 
trary to common expectation that ferromagnetic behav- 
ior is prop er of macroscopic systems only. Theoretical 
works [l2l . [isj , attempted to explain such ferromagnetic 
behavior in finite systems using microscopic models with 
on site anisotropy, inducing an effective energy barrier 



and thus leading to large magnetic reversal times and to 
ferromagnetic behavior due to finite measurement time. 

In this paper we take a different approach modeling 
both anisotropy and long-range with some suitable spin- 
spin interaction toy model as in [J but, differently from 
there, and in order to produce results closer to real ex- 
periments, we put the system in contact with a thermal 
bath. Despite its simplicity, it can be easily fitted to more 
physical models: for instance spin systems with dipole 
interaction in 3-D, have both long range and anisotropic 
spin-spin interactions. 

There are many different ways to model a thermal envi- 
ronment, especially when thermalization of a long range 
system is needed. Here we take the simplest route and 
assume that the environment is able to produce a Gibbs 
distribution for the system energies. 

In [3] we have found, for isolated anisotropic systems 
with an easy axis of magnetization (defined as the di- 
rection of the magnetization in the ground state energy 
configuration) that the constant energy surface is discon- 
nected in two regions with positive and negative magne- 
tization along the easy axis. This disconnection occurs 
below a critical energy threshold which has been called 
Topological Non-connectivity Threshold. Since the phase 
space is disconnected we have exact ergodicity breaking 
and no dynamical trajectory can visit the whole constant 
energy surface, but is, instead, limited to the region in 
which it started. Moreover, being defined for all finite iV, 
where N is the number of spins, the ergodicity breaking is 
not related to the thermodynamic limit. The Topological 
Non-connectivity Threshold is an example of topological 
singularity, well studied recently [ll], and its existence 
has been pointed out in Ising models too[l|. Also, an ex- 
perimental test of ergodicity breaking has been proposed 

Even if the connection between ergodicity breaking and 
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anisotropy has been shown to be a general one, indepen- 
dent of the range of interaction among the spins, long 
range interacting systems behave quite differently from 
short ones in the thermodynamic limit. This consequence 
has been studied in where it was found that the discon- 
nected portion of the energy spectrum remains finite in 
the thermodynamic limit for long-range interacting sys- 
tems only, while it becomes negligible, in the same limit, 
for short-range interacting ones. 

The plan of the paper is the following: In Sec. II, we 
review and extend the results obtained in Ref . 0] con- 
cerning the microcanonical behavior of long range inter- 
acting systems. In Sec. Ill we present a detailed calcula- 
tion of the density of states, using large deviation tech- 
niques, for a Mean Field Hamiltonian and we compared 
it for a generic long range interacting system. Finally in 
Sec. IV we study the magnetic reversal time in the canon- 
ical ensemble. In particular, we show how the magnetic 
reversal time for the open system can be obtained from 
that of the isolated one (a non trivial result). Also, a 
very simple approximation allows to interpret the Etnt 
as a real energy barrier for many body spin systems at 
sufficiently low temperature. 



II. THE TNT: MICROCANONICAL RESULTS 

Our paradigmatic anisotropic model, with an ad- 
justable interaction range, which presents ergodicity 
breaking for any N , is described by the following Hamil- 
tonian: 



H = J 



^ QX QX 



QZ QZ 



(1) 



where Sf,Sf,Sf are the spin components, assumed to 
vary continuously; i,j = 1, ...A^ label the spin positions 
on a suitable lattice of spatial dimension d, and is the 
inter-spin spatial separation. Each spin satisfies l^l = 1 
and < a < oo parametrizes the range of interactions: 
decreasing range for increasing a, so that a = corre- 
sponds to an all-to-all interacting model (close to phe- 
nomenological anisotropic models), while a = oo refers 
to a nearest neighbor interacting spin model. — 1 < ?/ < 1 
parametrizes the degree of anisotropy and for rj — —1 the 
Hamiltonian ([1]) docs not have a single easy axis. 

In Eq. IH) the constant J > has been added in order 
to fix the scale of time and to describe the model as 
ferromagnetic. 

Needless to say, this is not the most general spin Hamil- 
tonian giving rise to a TNT, even if it can be considered 
the simplest non integrable Hamiltonian with a suitable 
energy threshold above which the phase space is divided. 

The minimum energy configuration, with energy Emm^ 
is attained when all spins are aligned along the Z— axis ^ 
which defines implicitly the easy axis of magnetization. 



The phase space for E = Emin contains only two spin 
configurations, parallel or anti-parallel to the Z— axis. 
Therefore, the phase space at the minimal energy is dis- 
connected, due to the uniaxial anisotropy and it consists 
of two points only. We may ask now when and whether at 
a higher energy the constant energy surface is connected. 
To this purpose, let us define the TNT energy Emt as 
the minimum energy compatible with the constraint of 
zero magnetization along the easy axis of magnetization 
(hereafter we call m the magnetization along the easy 
axis) 



Et. 



Min 



H 



It is clear that if Etnt > Emm the phase space will be 
disconnected for all energies E < Etnt- We call this sit- 
uation Topological Non-connection, and, as will become 
clear in a moment, its physical (dynamical as well as 
statistical) consequences are rather interesting. Indeed, 
since below the TNT the phase space is disconnected, no 
energy conserving dynamics can bring the system from a 
configuration with to > to a configuration with to < 0, 
thus indicating an ergodicity breaking (impossibility to 
visit the whole energy surface). 

A useful quantity measuring how large the discon- 
nected energy region is, compared to the total energy 
range, can be introduced i4||: 



Etnt Eniin 
Eniin 



(2) 



In 3 it has also been shown that the disconnection ratio 
r, for N c». 





const ^ 



for 

for 



a>d 
a < d. 



(3) 



where d is the dimension of the embedding lattice. Since 
this point has not been remarked in Ref. , let us stress 
here that the existence of a phase transition for a < d can 
be inferred from the finiteness of r in the thermodynamic 
limit. Indeed, for long range systems, in order to define 
the thermodynamic limit it is convenient to make the 
energy extensive. This can be achieved by multiplying 
the Hamiltonian by N /\Eniin\- If we define the energy 
per particle, e = E/N, we can write: 



A^ 



e-tnt 



\E„ 



Et. 



En 



N 



N 



Since below etnt the most probable magnetization is for 
sure different from zero, then the specific energy at which 
the most probable magnetization is zero will be greater 
than the minimum energy in the thermodynamic limit, 
thus implying a phase transition. On the other hand, let 
us also notice that when r neither the existence nor 
the absence of a phase transition can be deduced. 
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An estimate for the TNT was also given for a < d and 
large N. More precisely, it can be shown that [J]: 



4^™„ - E.rmn for 



for 



■n > Vcr 



(4) 



where 



1 - 2"/'' < 0, 



and -B'j^j^ is the minimal energy for a system of N/2 
spins. For r/ > rj^r the TNT is given by the minimum of 
the second term in Eq. ([T]) under the constraint m ~ 0, 
while for rj < rjcr it is given by the minimum, under the 
same constraint, of the first term in Eq. ([T]). For rj < 
and finite N, there is a competition between the two 
different TNTs, therefore, in what follows, we will fix the 
anisotropy parameter rj = 1. Needless to say this choice 
does not affect the generality of our results. 

Due to the disconnection, below Etnt the dynamical 
time of magnetic reversal is infinite while above and 
close to the energy threshold (for chaotic systems), it 
was found to diverge as a power lawQ: 



{E - Et: 



~ e 



AS 



P 

J rr. 



Po 



(5) 



where AS' is the entropic barrier between the most likely 
magnetic states. Here, PE{m) is the probability distri- 
bution of the magnetization m at fixed energy i?, so that 

Pmax = Max„ PE{m), 

and Po = Psim = 0). 

The divergence found in [3| also shows that the phase 
space becomes connected at Etnt, a non trivial result, 
which cannot be deduced from the true existence of Etnt- 

Another important result found in Q is that, for all-to- 
all interacting spins (a = 0), the exponent 7 oc A^. This 
is related with the extensivity of the entropy S{E^ m) = 
InPE(m) (here we set fc^ = 1) and explain the huge 
metastability of such states even for small systems (say 
A^ ~ 100) and not necessarily close to the threshold Etnt- 

We numerically checked that, even for other power-law 
decreasing potentials in the long-range case a < d, 

• a power law divergence at Etnt still occurs as given 
byEq. ©; 

• reversal time is still proportional to Pmax/ Po (same 
equation) ; 

• the exponent 7 = A^. 

In order to do that we computed PE{fn = 0) and Pmax 
for different systems using the Wang-Landau algorithm 



In Fig. [l]we show the power law divergence of Pmax /Po 
for different a and A^ values. In order to improve the 
presentation we choose as a variable on the A— axis 



E-Et:, 



X ^ 



E.. 



Et. 



(6) 



where Egtat has been defined as the energy at which 
Pmax = Po (that is when the probability distribution 
of the magnetization has a single peak). That way all 
curves have a common origin. 
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FIG. 1: \ogPmax/Po vs \og X for different A'^ values as indi- 
cated in the legend and (a) a = 0.5; (b) a = 1. 

When the isolated system has a chaotic dynamics we 
computed the magnetic reversal time from the direct in- 
tegration of the equations of motion and we compare it 
with the "statistical" time Pmax/Po as given by Eq. ([5]). 
We show this comparison in Fig. [5] where each point on 
the graph has a X coordinate Pmax/Po and a Y coordi- 
nate r. The straight lines indicated proportionality over 
3 orders of magnitude. 

The linear dependence 7 = A^ found in [3] for the case 
a = holds for generic a too. 

In Fig. [3] we show the results of our numerical simu- 
lation for a = 0.5,0.9, 1 . Each point 7, at fixed A^ has 
been obtained computing the statistical reversal time for 
different energies, as plot in Fig. [U using the power law 
([5]). Assuming a power law dependence 7 oc A^'^ we have 
found a ~ 1 (within numerical errors) for all cases a <C 1 
(we show for simplicity only the case a — 0.5 in Fig. [3]). 

On the other hand, for a ^ I we have numerical evi- 
dence of a slower dependence on A^: 7 ~ N"' with a < 1. 
In the same Fig. [3] we show for sake of comparison the 
critical case a — d = 1 where a = 0.78(2), and the close- 
to-critical case a = 0.9 where we have found a = 0.87(2). 
Even if these results indicate that the simple linear re- 
lation 7 oc A^ is valid for long-range interacting systems 
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FIG. 2: Dynamical reversal time r vs the statistical one 
Pmax/Po for a = 1 and different A'^ values as indicated in 
the legend. Straight lines are r = {l/k)Pmax /Po with k = 10 
for A'^ = 16 and k = 2 for A*' — 48. They have been drawn with 
the only purpose to guide the eye showing the proportionality 
between the two quantities over 3 orders of magnitude. 




FIG. 3: log 7 as a function of logA'^ for different a values. 
Full circles stand for a = 0.5 and dashed line is the linear 
fitting with slope 1.02(1). Open circles is the critical case 
a = d = 1 and full line is the linear fit with slope 0.78(2). 
Open squares are for a = 0.9. Standard fit procedure gives 
a — 0.87 thus signaling the presence of a smooth transition 
at the point a = 1 for finite A'^. 



only, care should be used to extend the results of the case 
7 = d = 1 to large since finite A'' effects are huge in 
this case. Numerical evidence for cr < 1 has also been 
found in the short range case (a 1 ) but it will be 
discussed elsewhere. 



Let us also assume that are random variables uni- 
formly distributed in [—1, 1]. 

We can compute the entropy per particle as a function 
of rriz using Cramer's theorem [24| . so that we have: 



III. DENSITY OF STATES 

The density of states (DOS) for a Mean-Field approx- 
imated Hamiltonian can be computed analytically, using 
large deviations techniques [Hi. In particular wc will 
show that p{E) :^ {E — Emin)^ for E close to E,nin- We 
will also give numerical evidence that this law still consti- 
tutes an excellent approximation of the full Hamiltonian 
([T]) and for generic power law interaction a ^ 0. 

Let us consider the following Mean-Field Hamiltonian: 



mf 




(7) 



which can be considered a Mean-Field approximation of 
the Hamiltonian ([!]), for low energy and a = 0. Defining 



N 

k 

and e = Hmf/N, Eq. ([7|) can be rewritten as 
e = —Jnvt. 



s{mz) = -supx [Xruz - InV'(A)] 



(8) 



where 



2A 



Taking the sup in Eq. ^ we get: 



(9) 



which defines A as a function of to^. It is easy to see 
that for nriz ^ 1 then A — > oo (we could as well consider 
7712 —1 of course, and the result would be the same) 
so we restict our considerations to \mz\ — 1. Simplifying 
the expression of ij) we have 



- 2A' 



and inverting Eq. 



WLz 
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FIG. 4: The specific entropy vs energy, obtained numerically 
for A'^ = 24 and different q values (symbols indicated in the 
legend) is compared with that of the Mean Field Hamiltonian 
(full curve) and with the power law (dashed line), see Eq. 
In the inset the entropy is shown vs the energy for the case 
TV = 24, a = 0.5. 



Metropolis algorithm introduced in [23| ■ The idea behind 
is performing a random walk in phase space within an en- 
ergy range defined by the system temperature. The prob- 
ability -P(-E, T) of visiting a configuration with energy E 
and temperature T, obtained keeping a histogram of the 
energy values found during a Metropolis run, is related to 
the DOS through the Boltzman factor exp(--B/r). That 
provides us a conceptually simple way of determining the 
DOS. However, due to finite run time, P{E,T) will only 
contain information near (E) and we must combine the 
results for runs at different temperatures to obtain the 
complete DOS over the entire energy range. 

As one can see in Fig. H] the entropy per particle, in 
the long range case, is almost independent of the ran ge o f 
the interaction, also confirming a result obtained in [25| . 
Moreover, the theoretical approach gives a very good ap- 
proximation of the entropy per particle at low energy. 
When the energy is increased, the first term in the full 
Hamiltonian ([1]) becomes important and some deviations 
appear (see for instance the upper right corner in Fig. 
|4]). Needless to say the excellent agreement confirms the 
power law behavior for the DOS, Eq. pip , even for en- 
ergy values sufficiently high. As an example in the inset 
of Fig. m we can see that the power law expression (fTTj) 
holds up to Etnt- 



From Eq. ([8]) we obtain: 

s{6) —lnS + const. = In ( ^^^^ j -I- const., (10) 

\ Cmin / 

since 

e = ~Jml ~ - J(l - 2(5) = e™„ + 26 J, 

and Cmin — — J. From that wc immediately have that at 
low energy, 

piE) ~ (-B - Ermnf. (11) 

The next leading order in S can be calculated from 
e = -Jml = -J(l -<5)^ 

so that 
and 



s(e) = In 1 - J + const. (12) 

It is immediate to see that Eq. ([T^ gives Eq. (|10p for 
e ~ Cmin- We compared this analytical result for the 
Mean-Field model with the DOS computed numerically 
for the full Hamiltonian ([T]) and different a values in 
Fig. m The DOS has been calculated using a modified 



IV. MAGNETIC REVERSAL TIME 

Since the TNT has been introduced for isolated sys- 
tems, question arises if and how it can be defined when 
the system is in contact with a thermal bath. From the 
theoretical point of view we might expect that due to 
thermal noise the magnetization will be able, soon or 
later, to change its sign at any temperature T, thus sup- 
pressing the ergodicity breaking found in isolated sys- 
tems. Therefore, strictly speaking, a critical temperature 
below which the phase space is topologically disconnected 
for open finite systems does not exist. 

Nevertheless we are here interested in more practical 
questions, for instance: Will the energy threshold Emt 
still determine the magnetic reversal time in presence of 
temperature as it does in isolated chaotic systems? Can 
we predict the dependence of reversal time from temper- 
ature or any other system parameters, like the number 
of particles? 

Since the system is in contact with a thermal bath we 
may consider it as a member of a canonical ensemble. 
We may properly define the probability density to have 
a certain energy value E at the temperature T: P{E, T). 
Considering all members of the ensemble as independent 
objects we may guess that when the average energy {E) is 
much less than Emt and the probability density P{E, T) 
sufficiently peaked around its average value, the majority 
of the members of the ensemble will not cross the barrier, 
or at least, the probability of crossing it will be small. 
On the other hand for mean energy {E) on the order 
of Etnt each member will be allowed to jump, with a 
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time essentially given by the microcanonical expression 
Eq. dSD. 

Let us further assume, following the standard fluctu- 
ation theory [T^, [l^, that the magnetic reversal times 
between states with opposite magnetization are deter- 
mined by the free energy barrier between states at 
the most probable magnetization and states with zero 
magnetization: 



where 



fAA\ Max,„ Pt(to 

r cx exp — — — — . 

^ \ T J PT{m = 0) 



PtM = exp[-A(T,m)/T], 



(13) 



is the probability density to have magnetization m at the 
temperature T. Since MaxmiPT ("m)] is usually a slow 
varying function of the temperature, we can write 

1 



Pj.{m = 0) 



The crucial point now, is to obtain such value using the 
microcanonical results obtained in the previous Section, 
namely: 



Prim = 0) = ^ / PE{m = 0) e"^/^ p(E) dE, (14) 
Zt J 

where p{E) is the density of states and 

Zt = y e-^/^ piE) dE, (15) 



is the partition function. Since Pe{^ = 0) = for 
E < Etnt, the ergodicity breaking acts as a cut-off en- 
ergy of the integral p4|) , so that if the average energy is 
well below Etnt we can expect very large average reversal 
times. 

In order to verify that Eq. (fT4|) actually gives the mag- 
netic reversal time, we simulated the dynamics of a spin 
system in contact with a thermal bath in two differ- 
ent ways, the Metropolis algorithm [l^, and using the 
stochastic differential equations of the Langevin type as 
suggested in [2]| . 

In the Metropolis dynamics the change in the spin di- 
rection has been taken at each step completely random 
on the unit sphere while in the Langevin approach a small 
dissipation has been added. We checked that the results 
are independent of such dissipation and that the two ap- 
proaches give the same results. 

Obviously, both approaches can give directly the dis- 
tribution function: PT{m = 0), but we prefer here the 
direct calculation of the density of states and thus the 
possibility to get Pt(to = 0) for any temperature [23|, 
with less numerical effort and greater reliability. As for 
the reversal time it is quite obvious that an arbitrary 
multiplicative factor should be considered when the two 
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FIG. 5: Average reversal time r as a function of the rescaled 
temperature T/J for different interaction range, a) all-to-all 
a = 0, A = 24; b) long-range case a = 0.5, = 24; c) 
critical-case, q = 1, A'^ = 20. Circles are numerical data, 
dashed line is the integral calculated in Eq. (I14|l . 



different approaches are compared (the unit of time in 
the Metropolis approach is given by a random spin flip). 

Results are shown in Fig. O where the average rever- 
sal time (obtained with the Metropolis dynamics) vs the 
rescaled temperature T/J has been considered for differ- 
ent a and A^ values as indicated in the caption. As one 
can see the agreement between the integral ([T4|) (dashed 
line in Fig. [5]) and the numerical results (full circles) is 
excellent over 4 order of magnitude. It is worth of men- 
tion that no parameter fitting, other than a multiplicative 
constant has been used. 

It is also remarkable that a small variation in the tem- 
perature scale (about a factor 2) generates a huge vari- 
ation of the average time (roughly 3-4 orders of magni- 
tude). This signals an exponential dependence by the 
inverse temperature. Nevertheless a simple temperature 
dependence cannot be found in general, even if all curves 
in Fig.Ocan be fitted by the function T°- exp{b/T), with 
a and b fitting parameters. 

The exponential 1/T dependence can be understood 
in the limit of low temperature. Unfortunately we can- 
not compute directly the average time for temperatures 
lower than those shown, due to computer capability, even 
though we can study the asymptotic behavior of the in- 
tegral (|14p for low temperature. 



In order to obtain the low temperature behavior of 
the magnetic reversal times, we can use the saddle point 
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approximation to Eq. (|14|) getting: 

Prim = 0)^ PE'{m = 0) e(<^>-^*)/^ ^SiE')-Si(E)) ^ 



where {E) and E* are given by: 

1 dS ds 



where S = Inp, and 

T ^ dE 



(El 



N 



E* - Et. 



(16) 
(17) 

(18) 



Even if Eq. (|24p has been obtained for low temperature 
(T ^ Tcr), it should be kept in mind that this is a clas- 
sical model so that for T ^ 0, when quantum effects 
become important, it looses its validity [l^). 

The law has been checked numerically in Fig. 

|6] where the integral has been calculated for very 
low temperature and compared with the true exponen- 
tial law. As one can see asymptotically they are very 
close over many order of magnitude. 



where we used Pe* (m = 0) cx {E* - Etnt) , see Sect. II. 
An approximate expression for p?)) and (fT8|) can be ob- 
tained for small temperature. Indeed, using for the en- 
tropy the expression Eq. obtained in Sect. Ill, and 
inverting Eq. p7p . one obtains: 



(19) 




so that, for T <^ |emi„|/2, we get: 



< E>= E„ 



NT ( NT 

l + - +0 



Eri 



E„ 



(20) 



In the same way Eq. (jlSp can be written as: 



E* = Et: 



NT 



NT ^fNT\- 



(21) 



where A = Et^t — E^,^in and, for temperature sufficiently 
low, 



2N 



2 ' 



(22) 



we have 



E* = E. 



NT. 



Eq. (|16p can be further simplified, using the approxi- 
mated expressions for {E) and E* obtained above and 
Eq. lU) for the DOS: 



S{E*) 

sm) 



N ln{A + NT), 
N InTVT, 



(23) 



PB*(m-0) ~ {E*~Etni 



(NT) 



N 



Finally, neglecting terms of order {NT/ A): 



1 



Prim ^ 0) 



exp 



Etnt — E„ 

T 



(24) 
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FIG. 6: Average reversal time r calculated from Eq. (|14|l as a 
function of the rescaled temperature Tcr/T for different inter- 
action range a = (full circles) and a — 0.5 (open circles) and 
N — 24. Dashed and dotted lines represent their asymptotic 
value, as given by Eq. (|24p . 



Eq. ((24|) represents the central result of this paper. 
Indeed it allows to compute directly the reversal time in 
presence of a thermal bath at low temperature T without 
any complicated statistical calculations but the knowl- 
edge of the Hamiltonian itself. Moreover, the calculation 
of both the ground state energy and the Topological Non- 
connectivity Threshold constitutes a mechanical problem 
and they can be easily estimated even for complicated 
models. 

Furthermore, it also has some suggestive interpreta- 
tion. If we consider the path followed by the magnetiza- 
tion, as a random path of a Brownian particle between 
two potential wells separated by a potential barrier AU, 
according to Kramer's theory [13, [3 the average tran- 
sition time between the two wells follows the Arrhenius 
law: 



T - exp {AU/T) . 

Therefore it is clear that the disconnected energy region 
A = Etnt — Emin Can bc thought of as the real potential 
barrier felt by the magnetization. 
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In same way, the critical temperature Tcr has the physi- 
cal meaning of the specific energy barrier. It is interesting 
to note that the condition T <^ Tcr is not too restrictive, 
at least for long-range systems. Indeed, taking into ac- 
count that for large N [J]: 

lE^iJ 2-2" 1 

2N -2(2-a)(l-a) ' 

only at criticality (a = 1) it does not depend on N (and 
Tcr = In 2), while generally it grows with the number of 
particles. This is not at all surprising for a long range 
systems; indeed, if we make the energy of the system 
extensive, (multiplying the Hamiltonian by N /\E,-ain\)^ 
we have Tcr — ''/2, which is finite for any interaction 
range. 

Last, but not least, let us remark that the model given 
by Hamiltonian ([1]) at criticality {a = 1) is very interest- 
ing. Indeed, for a — 1 the parameter r (the ratio between 
the disconnected portion of energy space compared to the 
full one) goes to zero in the thermodynamical limit. The 
difference with the short range case is that it goes to 
zero logarithmically, instead of a power law: r ~ 1/ IniV. 
This simply means the existence of an effective phase 
transition for finite systems at criticality. 



V. CONCLUSIONS 



In conclusion we have shown that signatures of the 
topological disconnection persist when a long range in- 
teracting spin system is put in contact with a thermal 
bath. More precisely for temperature sufficiently low we 
recover the Arrhenius law for the magnetization rever- 
sal time T oc exp(A/T) similar to the reversal time for 
a Brownian particle jumping across a potential barrier 
A = Etnt — Emin- In other words the magnetization be- 
haves as a stochastic variable and the potential barrier 
is exactly given by the energy distance between Etnt and 
Emin- This proves the exponential dependence of the 
reversal times from the number of particles and stable 
ferromagnetism even for small systems with long range 
interaction and room temperature. The results presented 
in this paper can be experimentally verified, using for in- 
stance the physical system discussed in [15j, or in a 3-D 
spin system with dipole interactions. 
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